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Abstract 

The adenosine A3 receptor (AAR) belongs to a subfamily of G protein-coupled receptors and 
is an important therapeutic target for conditions including inflammation and cancer. The clinical 
compounds CF101 and CF102 are potent and selective A3AR agonists, but the structural basis 
of their recognition was unknown. Here we present the cryogenic electron microscopy 
structures of the full-length human AsAR bound to CF101 and CF102 at 3.3-3.2 A resolution 
in complex with heterotrimeric Gi protein. These agonists bind within the orthosteric pocket, 
with their adenine components engaging in conserved interactions while their substituted 3- 
iodobenzyl groups exhibit different orientations. Swapping extracellular loop 3 (ECL3) of AsAR 
onto other adenosine receptor subtypes enabled CF101/CF102 binding and receptor activaton, 
and mutations in key residues, including His??7, Ser*^? and Ser$*?? that form a unique 
subpocket in A3AR, abolished receptor activation, highlighting these structural elements are 
critical for ligand selectivity. Compared to inactive AzaAR, the AsAR structures reveal 
conserved mechanism of receptor activation, including an outward shift of TM6. These 
structures provide key insights into molecular recognition and signaling mechanisms of A3AR, 
which should aid rational design of subtype-selective ligands targeting this important class of 
adenosine receptors. 
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Introduction 

The adenosine receptor subfamily of G protein-coupled receptors consists of four subtypes: 
A1, Aza, A28, and Asl: 2]. These receptors are activated by the endogenous ligand, adenosine, 
to transduce downstream signals that mediate a number of important physiological and 
pathalogical roles including immunomodulation, energy balance, cardiac function, 
neuroprotection, etc 5). AsAR is expressed in various tissues including the brain, heart, lungs, 
liver, kidneys, and immune cells!!. Through its signaling functions, AsAR participates in 
regulating cardiac function, vasodilation, inhibition of inflammation, protection against 
ischemia-reperfusion injury, and suppression of oxidative stress. Additionally, AsAR is highly 
expressed in a number of tumor cells, suggesting its potential as a therapeutic target for 
suppressing cancer cell proliferation!®®). 


The activation of Az4AR and AæAR predominantly elicits stimulatory G protein (Gs) signaling, 
while AiAR and AsAR exhibit a preference for coupling to inhibitory G protein (Gi), leading to 
inhibition of adenylate cyclase and decreased intracellular cyclic AMP!!. Based on the 
chemical structure of adenosine, numerous agonists and antagonists have been designed and 
tested against AAR for disease indications such as cancer, inflammation, and pain?! . 
Previous studies suggest that modifications at the N position on the purine and 5'-N position 
on ribose group based on the adenosine framework yield potent A3AR agonists with high 
subtype selectivity! ®-121. CF101 and CF102 are representatives of such modification strategy 
with similar nucleoside core structure and only one chloro-substituent difference, both 
demonstrate high affinity and selectivity for ASARÜ*5, CF101 showed efficacy in Phase Ill 
trials for psoriasis and rheumatoid arthritis? while CF102 is in clinical trials for hepatocellular 
carcinoma and non-alcoholic steatohepatitis (NASH)!"® 171, 


Adenosine receptors mediate many important functions and their wide expression makes 
subtype selectivity of ligands more critical to reduce possible side effects !"* 191. Elucidating the 
structural basis of selective ligand binding to AsAR could help in improving the drug design 
process. To date, no AsAR structure has been reported. In this paper, we present the cryo-EM 
structures of A;sAR bound to Gi in the presence of CF101 and CF102, which reveals the basis 
of ligand recognition and ligand-induced activation mechanism of AsAR. Our works provide 
important insights for designing effective A3AR-targeted therapies, and more broadly for the 
subfamily of adenosine receptors. 


Results and Discussion 

Overall structures of the complexes 

CF101 and CF102 are A3AR agonists containing modifications to the ribose and adenine 
moieties that confer potent binding to AsAR (Fig. 1a). Specifically, CF101 and CF102 have a 
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5'-N-methylcarboxamide substitution on the ribose group and a NS-(3-iodobenzyl) substitution 
on the adenine base (Fig. 1a). These modifications result in significantly higher ASAR potency 
compared to the endogenous As3AR agonist adenosine. We verified the selectivity of these 
nucleoside-derived compounds for A3AR versus other adenosine receptor subtypes 
(Ai/Aza/Azp-AR) using NanoBiT association assays (Fig. 1b-d). While adenosine activated four 
subtypes with similar micromolar potencies, CF101 and CF102 displayed strongest potency 
of ^3 nM on A3AR but had weak or negligible response on other subtypes of adenosine 
receptors. 


We used NanoBiT tether strategy to stabilize the full-length ASAR-G protein complexes, as it 
has been used for many GPCR structural studies!*°-2] (Supplementary Fig. 1). AsAR used in 
this study had an N-terminal thermostabilized apocytochrome b562RIL (BRIL) fusion to 
enhance its expression, which is co-expressed with G protein subunits and scFv16, an 
antibody fragment that is used to further stabilize the receptor G protein complex. For the 
CF101-A3AR-Gi complex, data from 20,779 movies comprising 271,323 particles was used to 
determined the structure at 3.29 À resolution (Supplementary Fig. 2, Supplementary Table 1). 
For the CF102-AsAR-Gi complex, data from 13,581 movies yielding 283,561 particles was 
used to determined the structure at a resolution of 3.19 À (Supplementary Fig. 3, 
Supplementary Table 1). The structures of the CF101/CF102-AsAR-Gi complexes revealed 
that the ligands occupy the orthosteric binding pocket, with the core structures modeled clearly 
into the cryo-EM density at the center of the receptor transmembrane helicals (TMs) (Fig. 1e- 
h). 


The structures showed the canonical seven-transmembrane architecture for A3AR, with the 
intracellular domains occupied by the a5 helix of Gai for Gi coupling. The density maps enabled 
modelling of most of the structures, except for AsAR N-terminus residues M1-L8, third 
intracellular loop N21 1-Y222, C-terminus V301-E318, and the alpha-helical domain of Gai. The 
extracellular loop M151-S165 was also less defined but the backbone could be established 
(Supplementary Fig. 4). Aside from these regions, the models were well-resolved. Overall, the 
two agonist-bound complexes were highly similar, with 0.593 À root mean square deviation 
(RMSD) for the whole receptor. 


Binding mode of CF101/CF102 in A3AR orthosteric site 

The AsAR agonists CF101 and CF102 bind at conserved orthosteric pocket formed by ECL2, 
TM3, TM5, TM6 and TV7, akin to the endogenous ligand adenosine bound to other adenosine 
receptor subtypes (Fig. 2a-b). However, the orientations of the modified 3-iodobenzyl moieties 
differ between CF101 and CF102. The adenine core mediates conserved receptor interactions 


commonly seen in other adenosine receptors": 23. 24, Notably, the adenine pyrimidine forms 
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Tr-stacks against F^55?, and the F^95?A mutation greatly affected the ability of CF101/CF102 to 
induce the receptor activation (Fig. 2c-f, Supplementary Table 2). Additionally, ribose and 3- 
iodophenyl groups form hydrogen bonds with polar side chains at positions 3.36, 6.55 and 
7.43, which are key for recognition of nucleoside ligands by all adenosine receptors (Fig. 2c-f, 
Supplementary Table 2). 


The ligand binding pocket is mainly comprosed of hydrophobic residues, including position 
3.33, 5.38, 5.47, 6.48, 6.51 and 7.39, which form hydrophobic contacts that are important for 
CF101/CF102 potencies (Fig. 2c-f, Supplementary Table 2). Alanine mutations at these 
positions severely reduced agonists’ ability to induce receptor activation. His*?7 and Ser?^? 
participate van der Waals contacts with the bound ligands, their alanine mutations also 
affected activity (Fig. 2c-f, Supplementary Table 2). The side chains from M174535 and L2647?5 
in the receptor form hydrophobic interactions with the 3-iodophenyl group extended from the 
Nê position of the adenosine base of CF101. In contrast, the corresponding group of CF102 is 
surrounded by V1695°!? and L2647°5 from the receptor. Alanine mutations on these residues 
did not significantly affect the potency of the compounds on A3AR (Supplementary Fig. 6a, 
Supplementary Table 2), suggesting that the 3-iodophenyl substituents may exist alternative 
states at the receptor extracellular domains. This demonstrates that the NÓ position may 
accommodate various substituted groups through distinct conformations in the A3AR pocket. 


Moreover, CF102 is a 2-chloro derivative of CF101 (Fig. 1a). Y15135 in CF102-bound AAR 
forms hydrophobic contact with the 2-chloro group in CF102. Meanwhile, Y15135 in AsAR 
formed T-T interaction with Y2657?6 in TM7. The Y15'°5A mutation in AsAR abolishes the 
agonist activity of both CF102 and CF101 (Fig. 2c-f, Supplementary Table 2). According to the 
reports, modifications to the 2-position of the adenosine structure tend to well-tolerated in 
binding AsARU^I, whether it a small or large group, even the macrocycle group linked with the 
group from Nê moeity*!. Elucidation of the subtle structural variations in ligand and receptor 
interactions thus provides molecular insight into the conformational adaptability and binding 
pose governing molecular recognition at A3AR. 


ECL3 in adenosine receptors 

CF101 and CF102 show high selectivity on A3AR rather than other subtypes. Analysis the 
sequence of adenosine receptors reveals strong conservation within TMs, while the 
extracellular loops diverge among subtypes (Supplementary Fig. 5). ECL1 shows relatively 
distant from the orthosteric site. F1684^55? in ECL2 provides the key t-m interactions with each 
agonists in adenosine receptors. However, A3AR possesses a shorter ECL3 than other 
subtypes (Fig. 3a). The shorter ECL3 may rigidity ASAR to minimize its conformational changes 
for specific ligand binding. 
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To assess the role of ECL3 in A3AR, we engineered chimeric receptors by grafting ECL3 from 
A3AR onto the backbones of other adenosine receptors. The chimeric receptors gained the 
ability to bind CF101 and CF102 with increased efficacy or potency (Fig. 3b-c, Supplementary 
Table 3). These findings suggest that ECL3 could serve as a structural factor mediating the 
selective recognition CF101 and CF102 by AsAR. According to the reported structure-activity 
relationship of ligands at the AsAR, numerous N$-substituted adenosine derivatives were 
synthesized. Too small or bulky groups at the N® position would reduce its potency or affinity 
on AsARI?: 14, The N° position on adenosine is projecting outwards into the binding pocket of 
A3AR and is in close spatial proximity to the ECL3 region of the AAR. So the ECL3 loop is an 
important consideration in structure-activity studies of N$-modified adenosine derivatives 
against adenosine receptors. Delineating these subtle structural variations provides molecular 
insight into the selectivity of structurally analogous ligands for adenosine receptors. 


Binding pocket residues across adenosine receptors 

Furthermore, A3AR shared lowest identity with other subtypes among adenosine receptors. 
The sequence analysis reveals that the A3AR may confer selectivity through different residue 
types in the orthosteric binding pocket (Fig. 4a), which include residues at positions 3.32 
(L/V/V/V, the residue in As/Ai/AzA/Aeg-AR), 3.37 (H/Q/Q/Q), 5.42 (S/N/N/N), 5.47 (V/V/V), 
6.52 (S/H/H/H) and 6.58 (I/T/T/T) (Fig. 4b, Supplementary Fig. 7). These residues in AsAR 
were mutated to the corresponding residues of other receptor subtypes, with aim to assess 
the impact of these AsAR mutants on the activity of CF101 and CF102 (Fig. 4c-d). 


The leucine at position 3.32 in AsAR, versus valine in other subtypes, did not affect 
CF101/CF102 activity when mutated to valine, consistent with similar hydrophobic properties 
(Fig. 4c-d, Supplementary Fig. 7, Supplementary Table 2). Likewise, mutating isoleucine at 
positions 5.47 and 6.58 in AAR to the valine and threonine found in other subtypes slightly 
impacted activation (Fig. 4c-d, Supplementary Table 2). This suggests that the slightly shorter 
side chains present in valine and threonine do not impair binding. 


The side chain of H937 and S5^? in AsAR form a hydrogen bond , which could not be formed by 
corresponding residues Q??7 and N54? in other adenosine receptor subtypes (Fig. 4e-i). The 
H?3'Q mutation showed a limited impact on the activity of CF101/CF102 but the S°4*N 
mutation almost abolished the ability of CF101/CF102 to induce the receptor activation (Fig. 
4c-d, Supplementary Table 2). Additionally, mutating S95? to histidine (H) also severely 
decreased CF101 and CF102 activity, likely due to unfavorable steric or electronic properties 
of the longer histidine side chain (Fig. 4c-d, Supplementary Fig. 7, Supplementary Table 2). In 
contrast to other adenosine receptor subtypes, H??/, S54? and S95? form a distinctive 


6 


189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
201 
208 
209 
210 
211 
212 
213 
214 
215 
216 
217 
218 
219 
220 
221 
222 
223 
224 
225 
226 


subpocket in AsAR to accommodate the 5'-N-methylcarboxamide from the ribose (Fig. 4e-i, 
Supplementary Fig. 8). These results implicate this microdomain serves as a structural 
determinant for stabilizing CF101 and CF102 in A3AR versus other subtypes. Together, these 
findings demonstrate that minor sequence changes in the receptors could impact their 
conformations, thereby affecting ligand binding specificity. 


The activation mechanism of adenosine receptor A3AR 

Structural comparisons between active, agonist-bound A3AR complexes and an inactive, 
antagonist-bound AzaAR structure (PDB ID: 4EIY)P?9 reveal hallmarks of conformational 
changes associated with GPCR activation” 281, The AsAR structures exhibit an outward 
movement of TM6 compared to inactive AzAAR, shifting 11.6 A based on measurements of 
residue Glu$39? at Ca atoms in receptors, analogous to movements seen in other activated 
class A GPCRs upon G protein coupling (Fig. 5a). Additional rearrangements of activation 
include inward movements of TM1 and TM7 and an upward shift of TM3 in AsAR relative to 
inactive AeAAR (Fig. 5b-d). The AasAR agonists CF101 and CF102 dock deeper into the 
orthosteric pocket compared to the AzaAR antagonist ZM-241385, enabling engagement of 
conserved GPCR activation motifs like the “toggle switch" on W848 and transmission switch 
motifs D949R950Y35! and N749P750xxY753 (Fig. 5e-h). These microswitches trigger the 
conformational changes from the ligand binding pocket to the cytoplasmic G protein coupling 
interface. The series of structural transitions led to the rearrangements of receptor to enable 
G protein coupling and activation. By providing near-atomic resolution of A3AR activation 
mechanisms, these findings reveal fundamental insights into the relationship between ligand 


recognition in receptor, and activation of downstream signaling. 


G protein coupling of adenosine receptors 

All adenosine receptors exhibit differential G protein coupling preferences that correlate with 
distinct conformational orientations of the associated G proteins': °° 24, A,AR and AsAR 
preferentially couple to inhibitory Gi proteins, while Az4AR and AesAR primarily couple to 
stimulatory Gs proteins to mediate intracellular signaling cascades (Fig. 6a). Structural 
comparisons reveal conformational differences between Gi- and Gs;-coupled adenosine 
receptors complexes. Specifically, the TM6 helix of AiAR/AsAR-Gi complexes shows a 3.1 A 
inward shift to accommodate Gi binding compared to A24AR/A2BAR-Gs complexes (Fig. 6b). 
Additionally, the a5 helix of Gas subunits displays a 8.6 À displacement relative to its 
orientation in Gi complexes based on measurements of the Ca atom of Ga" 59$ (Fig. 6c). The 
aN helix of Ga; exhibits a 3.3 A tilt compared to Gs when measuring the Ca of Ga" 38 (Fig. 6d). 
Overall, the structural arrangements of AiAR/A3AR-Gi complexes closely resemble each other, 
similar to the consistency observed between AzaAR/AzsAR-G; complexes. These findings 
reveal that Gi-coupled adenosine receptors adopt conserved G protein-binding conformations 
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that differ distinctly from those of Gs-coupled adenosine receptors. Elucidation of the structural 
differences governing adenosine receptors-G protein interactions provides fundamental 
insights into the molecular determinants of G protein coupling specificity. 


Conclusion 

In conclusion, we have determined the cryo-EM structures of the A3AR bound to selective 
agonists CF101 and CF102 with heterotrimeric Gi protein. Despite the conserved binding of 
the core adenosine moiety, the structures revealed differences in the orientations of the N° 
substituted groups in CF101 and CF102. We identified ECL3 and key pocket residues His? ?", 
Ser*^ and Ser®5? that confer selectivity over other adenosine receptor subtypes by 
mutagenesis studies. Comparison to an inactive A2AAR structure provided insight into the 
conformational changes associated with A3AR activation and G protein coupling. By 
elucidating the molecular mechanisms governing ligand recognition, signaling, and subtype 
selectivity, A3AR structures significantly advance our fundamental understanding of this 
important drug target. The findings pave the way for structure-guided design of improved 
therapeutics targeting adenosine receptors for the treatment of cancer, inflammation, and 
other diseases. 
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350 

351 Fig 1. Cryo-EM structures of CF101-A3AR-G; and CF102-A3AR-G; complexes. 

352 a. Chemical structures of the adenosine, CF101 and CF102. Highlighting modifications at the 
353 5'-N-methylcarboxamide in the ribose group and N® and C2 positions of the adenoisine group. 
354 The atom numbering was indicated in blue. CF101, is also named IB-MECA and N$-(3- 
355 iodobenzyl)adenosine-5'-N-methyluronamide. CF102, is also named CI-IB-MECA and 2- 
356 chloro-N®-(3-iodobenzyl)adenosine-5'-N-methyluronamide. b-d. NanoBiT association assays 
351 monitoring ligand activity on adenosine receptors for adenosine (b), CF101 (c) and CF102 (d) 
358 from three independent experiments in triplicate (n23). e, f. Cryo-EM map (e) and model (f) of 
359 the CF101-A3AR-Gi complex, with inset showing CF101 density. The density map in the inset 
360 is shown at 0.232 threshold. g, h. Cryo-EM map (g) and model (h) of the CF102-A3AR-G; 
361 complex, with inset showing CF102 density. The density map in the inset is shown at 0.17 
362 threshold. Subunits are colored as indicated. 
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Fig 2. The orthosteric binding pocket. 
a, b Detailed interactions between AsAR and CF101 (a) or CF102 (b) from the membrane 
plane. Residues involved in ligand interaction are colored blue and pink in two complexes, 
respectively. Black dashed lines indicate hydrogen bonds. c-f Dose-response curves of 
muntants of AsAR induced by CF101 (upper panels, c, e) or CF102 (lower panels, d, f) using 
NanoBiT assay. 
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Fig 3. Swapping ECL3 increases CF101/CF102 potency at adenosine A1/AzA/Ace receptor 
subtypes. 

a Superposition of ECL3 in the adenosine receptors shows AsAR reveals A3AR has a shorter 
ECL3. Other TMs were omitted. b, c Effects of CF101 and CF102 were tested on the A1AR, 
AzaAR, and AæAR and their corresponding mutants containing the swapped ECL3 from the 
A3AR using NanoBiT assays. 
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382 (7LD4) (6GDG) (8HDP) 
383 Fig 4. Key residues in the A3AR binding pocket. 
384 a Sequence alignment of the residues in the orthosteric binding pocket among the adenosine 
385 receptors. The conserved residues were colored with blue and stars marked. The unique 
386 residues in A3AR, different from other adenosine receptors subtypes were colored in orange, 
387 while the others subtpyes at corresponding position were colored in green. All residues were 
388 annotated based on GPCR Ballesteros-Weinstein numbering scheme. b Superposition of 
389 adenosine receptors, the unique residues in A3AR from the other adenosine receptors are 
390 shown as yellow ball. c-d Effects of CF101/CF102 on AsAR mutants containing swapped 
391 residues from other adenosine receptors by NanoBiT assay. e-i The binding cavities of the 
302 adenosine receptors are generated in PyMOL, the cavities were depicted as gray. In AsAR, 
393 His?37,Ser*^? and Ser®5 form a subpocket, whereas these positions are conserved as 
394 GIn??7, Asn*^? and His9*? in the other adenosine receptor subtypes, respectively (His, H; 
295 Ser, S; Gin, Q; Asn, N). Dash lines in h and i showed the hydrogen bonds between His??" 


396 and Ser®-42, 
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A,,AR-ZM241385(inactive) 
A,AR-CF 101 (active) 
A,AR-CF 102(active) 


toggle switch PIF motif DRY motif NPxxY motif 
Fig 5. A;AR activation mechanism. 


a-b Superposition of active ASAR-CF101/CF102 complexes (blue/pink) with inactive AzAAR- 
ZM241385 complex (grey, PDB ID 4EIY). c,d Comparison of extracellular (c) and 
cytoplasmic (d) views of active AsAR and inactive AzAAR. e-g Conformational changes in 
conserved motifs, including toggle switch, PIF, DRY and NPxxY, upon CF101/CF102 binding 


to A3AR relative to inactive AzaAR-ZM241385. Arrows indicate movement directions. 
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Figure 6. G protein coupling of adenosine receptors. 


a Comparison of G protein conformations in Ai/A3AR-Gi and AeA/AesAR-Gs complexes. b 


Conformational comparison of TM6 in adenosine receptors, referencing toggle switch W$^ in 


TM6 of receptor. c-d Conformational comparison of a5 helix and aN helix in G protein 


adenosine receptors-G protein complexes. Arrows indicate movement directions. 
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Methods 

Construct design 

The full-length gene coding human A3AR was synthesized (Synbio) and subcloned into 
pFastBac vector using CloneExpress II one step cloning kit (Vazyme Biotech). A hemagglutinin 
signal peptide and thermostabilized apocytochrome b562RIL (BRIL) were fused at the N- 
terminal of AsAR to enhance receptor expression. To enhance complex stability, a NanoBiT 
tethering approach was used where an LgBiT domain was fused to the C-terminal of the 
receptor!’ A dual maltose-binding protein was linked after LgBiT through a tobacco etch virus 
protease site (TEV site) for further cleavage. A dominant-negative mutant of bovine Gai 
containing G203A/A326SP?l mutations was generated to stabilize the heterotrimeric GaiBy 
protein. Rat GB1 was fused with a HiBiT at C-terminal for structural complementation of LgBiT 
to form a NanoBiT. The single-chain variable fragement scFv16 was applied to bind the Gay 
protein for stabilization. Gai, GG1-HiBiT, Gy, and scFv16, were cloned into pFastBac vector 
(Supplementary Fig. 1a), respectively. 


Protein expression and purification 

The recombinant A3AR, Gai, GB1-HiBiT, Gy, and scFv16 were co-expresed in Trichoplusia ni 
High Five insect cells using the Bac-to-Bac baculovirus expression system. High Five cells 
was co-infected with the baculovirus at a cell density of 3.5 x 108 cells per milliliter. Fourty- 
eight hours later, the infected cells were harvested and stored at -80 °C until used. 


For the purification of the CF101-A3AR-Gi complex, cells pellets was thawed and resuspended 
in Buffer A (100 mM NaCl, 20 mM HEPES, pH 7.5) supplemented with protease inhibitor 
cocktail (TargetMol). Cells was lysed by dounce homogenization followed by centrifugation to 
remove unsoluble materials. The pellets were respended in Buffer B (100 mM NaCl, 10 %(v/v) 
glycerol, 20 mM HEPES, pH7.5) supplemented with 10 mM MgCl», 5 mM CaCle, 0.2 mM Tris- 
(2-carboxyethyl)phosphine (TCEP, Hampton Research) and protease inhibitor cocktail. The 
complex was formed by rotating at room temperature for 1 h after addition of 25 mU/mL 
apyrase and 10 uM CF101. After incubation, the sample was solubliezed in 0.5 %(w/v) lauryl 
maltose neopentylglycol (LMNG, Anatrace) and 0.1%(w/v) cholesteryl hemisucinate (CHS, 
Anatrace) for 3 h at 4 °C. The supernatant was clarified by centrifugation at 100,000x g for 40 
min. The supernatant was incubated with dextrin beads 6FF (Smart-Lifescience) for 3 h at 4 °C. 
The beads were loaded onto a gravity column and washed with 20 column volumes of Buffer 
C (100 mM NaCl, 2 mM MgCle, 10 uM CF 101, 0.2 mM TCEP, 0.01 %(w/v) LMNG, 0.002 %(w/v) 
CHS, 20 mM HEPES, pH 7.5). The complex was eluted with Buffer C supplemented with 10 
mM maltose and further concentrated using 100 kDa molecular weight cut-off concentrator. 
TEV protease was added to the concentrated protein at 4 °C overnight to cleave dual maltose 
binding protein from fusion protein. After digestion, sample was loaded onto Superdex 200 
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Increase 10/300 GL column (GE Healthcare) with Buffer D (100 mM NaCl, 2 mM MgCle, 10 
uM CF101, 0.1 mM TCEP, 0.00075 %(w/v) LMNG, 0.00025 %(w/v) glyco-diosgenin, 
0.0002 %(w/v) CHS, 20 mM HEPES, pH 7.5). The desired fractions were pooled and 
concentrated to 5-8 mg/mL for cryo-EM sample preparation. 


Cryo-EM data collection 

Cryo-EM grids were prepared with the Vitrobot Mark IV plunger (FEI) set to 8 °C and 100% 
humidity. Three-microliters of the CF101-A3AR-G; complex were applied to glow- discharged 
Quantifoil R1.2/1.3 holey carbon grids. The sample was incubated for 10 s on the grids before 
blotting for 3.5 s (double-sided, blot force 1) and flash-frozen in liquid ethane immediately. The 
same conditions were used for the CF101-A3AR-Gi complex sample. 


For CF101-As3AR-Gi complex, three datasets comprising 20,779 movies were collected on a 
Titan Krios equipped with a Gatan K3 direct electron detection device at 300 kV with a 
magnification of 105,000 corresponding to a pixel size 0.824 À. Image acquisition was 
performed with EPU Software (FEI Eindhoven, Netherlands). We collected a total of 36 frames 
accumulating to a total dose of 50 e A? over 2.5 s exposure. 


For CF102-As;AR-Gi complex dataset, two datasets totaling 13,581 movies were collected on 
a Titan Krios equipped with a Gatan K3 detector at 300 kV with a magnification of 105,000 and 
pixel size of 0.824 A, using EPU Software (FEI Eindhoven, Netherlands). Thirty-six frames 
were collected over a 2.5-s exposure to a dose of 50 e À?. 


Image processing 

MotionCor2 was used to perform the frame-based motion-correction algorithm to generate 
drift-corrected micrograph for further processing, and CTFFIND4 provided estimation of 
contrast transfer function (CTF) parameters! 321, 


For the CF101-A3AR-G; dataset, the previously resolved structure of BAY 60-6583-A2BAR-Gs 
[21] was used as a reference for automatic particle picking in Relion 3.0033. Particle picking and 
extraction yielded 4,550,294 particles after 2D classification clearance, which were imported 
into cryoSPARCP!, Four rounds of 2D classification selected 1,267,837 particles, followed by 
two rounds of 3D heterogenous refinement giving 982,833 particles. After two additional 
rounds of 2D classification and two rounds of heterogenous refinement, 271,323 particles were 
refined to a structure at 3.29 À global resolution using non-uniform refinement (Supplementary 
Fig. 2). 
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For CF102-AaAR-Gi complex dataset, the BAY 60-6583-A2BAR-Gs structurel?!! was again used 
for reference-based particle picking. 4,090,959 and 4,833,382 particles were autopicked and 
extracted from Dataseti and Dateset2, respectively. For Dataseti, two rounds of 2D 
classification were used to separate out 1,070,085 particles. Masked 3D classification on the 
receptor part was used to separate out 175,747 particles that resulted to a clearer density of 
A3AR. For Dataset2, two rounds of 2D classification and two rounds of 3D classification were 
performed to separate out 246,392 particles. After clearance, the remained particles from two 
datasets were combined and subjected to alignment-free 3D classification. 283,561 particles 
were remained and transferred in cryoSPARC'1. One round of heterogenous refinement 
yielded a final 102,581 particle were refined to a structure at 3.19 À global resolution using 
non-uniform refinement (Supplementary Fig. 3). 


Model building 

An A3AR structure predicted by AlphaFold2 was used as the starting reference models for 
receptors building. Structures of Gai, GB, Gy and the scFv16 were derived from PDB entry 
7EZHP6! were rigid body fit into the density. All models were fitted into the EM density map 
using UCSF Chimera!” followed by iterative rounds of manual adjustment and automated 
rebuilding in COOT®8! and PHENIXP?I, respectively. The model was finalized by rebuilding in 
ISOLDEU! followed by refinement in PHENIX with torsion-angle restraints to the input model. 
The final model statistics were validated using Comprehensive validation (cryo-EM) in PHENIX 
and provided in the supplementary information, Supplementary Table 1. All structural figures 
were prepared using Chimera'7l, Chimera X'4"], and PyMOL (Schrödinger, LLC.). 


NanoBiT assay 

To monitor G protein interaction with AiAR, AzaAR, AeeAR or AAR upon agonist stimulation, 
a NanoLuc-based NanoBiT enzyme complementation assay was used (Promega). The C- 
terminus of AiAR, AzaAR or AAR was fused to SmBiT, while LgBiT was fused to the N- 
terminus of G proteins. The C terminus of A3AR was fused with LgBiT, and the SmBiT element 
was fused to the N terminus of G proteins. HEK293 cells were seeded at 4 x 10^ cells/well on 
96-well plates and co-transfected with AR-SmBit and LgBiT-G protein plasmids at a 1:1 mass 
ratio. After 24 hours, cells were replaced with 40 uL fresh culture medium without fetal bovine 
serum. Ten microliter Nano-Glo Live Cell reagent was added followed the manufacturer's 
protocol (Promega, N2011), and incubated at 37 °C, 5 % CO» for 5 min. Another 25 uL culture 
medium containing various concentrations of compounds were added and incubated at room 
temperature for 10 minutes before measuring bioluminescence using EnVision multiplate 


reader (PerkinElmer). 
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Cell-surface expression assay 

Wild type AiAR, AzaAR, A2BAR or A3AR gene was subcloned into pcDNA3.0 vector with an N- 
terminal human influenza hemagglutinin tag (HA-tag). HEK293 were seeded at 4 x 10^ 
cells/well on 96-well plates and transfected with wild type (WT) or adenosine receptor mutants. 
After 24 hours, cells were washed with PBS buffer, fixed with 4 %(w/v) paraformaldehyde for 
15 min, and blocked with 2 %(w/v) bovin serum albumin (BSA) for 1 h. Next, cells were 
incubated with the polyclonal anti-HA (Sigma, H6908) overnight at 4 °C, followed by 1 h with 
horseradish peroxidase (HRP)-conjugated anti-rabbit antibody (Cell Signaling, 7074S) at room 
temperature. After washing, 50 uL tetramethylbenzidine (Sigma, T0440) was added for 30 min 
before stopping the reaction with 25 uL TMB substrate stop solution (Beyotime, P0215). 
Absorbance at 450 nm was measured on a FlexStation Ill microplate reader (Molecular 
Devices). 


Statistical analysis 

All functional study data were analyzed in Prism 8 (GraphPad) and presented as means + 
S.E.M. from at least three independent experiments. Concentration-response curves were 
evaluated with a three-parameter logistic equation. ECso values were calculated using the 
sigmoid three-parameter equation. Significance was determined by one-way ANOVA followed 
by multiple comparisons test, and *P « 0.05 vs. wild-type (WT) was considered statistically 
significant. 
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Supplementary Fig. 1 The expression and purification of ASAR-Gi complex. 
a. Schematic diagrams of the expression constructs of A3AR and G81 using the NanoBiT 
tethering approach. AsAR and G81 fused with LgBiT and HiBiT, respectively. b. Size- 
exclusion chromatography profile of the CF101-A3AR-Gi complex. c. SDS-PAGE of the arrow 
indicated peak fraction in (b). d. Size-exclusion chromatography profile of the CF102-A3AR- 
Gi complex. e. SDS-PAGE of the arrow indicated peak fraction in (d). 
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Dataset 1 Dataset 2 Dataset 3 


Micrographs (7,229 movies) Micrographs (6,528 movies) Micrographs (7,022 movies) 
Motion correction Motion correction Motion correction 
CTF estimation CTF estimation CTF estimation 
Autopick Autopick Autopick 
(4,767,541 particles) (1,646,832 particles) (3,555,824 particles) 
2 rounds 2D classification 2 rounds 2D classification 
(1,474,931 particles) (1,428,531 particles) 


Transfer to Cryosparc [4.550.204 particles) 
4 rounds 2D classification 
(1,267,837 particles) 


2 rounds heterogeneous Refinement 
(982,833 particles) 


2 rounds 2D classification 
(750,750 particles) 


2 rounds heterogeneous Refinement 
(430,696 particles) 


271,323 P 


Non-uniform Refinement 
d e 
(3.29 A) 


GSFSC Resolution: 3.29A 


Supplementary Fig. 2 Cryo-EM data processing of CF101-A3AR-G; complex. 

a. Representative image from cryo-EM dataset. Scale bar, 50 nm. b. Respresentative 2D 
average classification classes. Scale bar, 10 nm. c. Flow-chart of the cryo-EM data 
processing. d. FSC curves. e. The local resolution map. 
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Dataset 1 Dataset 2 
Micrographs (5,720 movies) Micrographs (7,861 movies) 


Motion correction Motion correction 
CTF estimation CTF estimation 
Autopick Autopick 
(4,090,959 particles) (4,833,382 particles) 

2 rounds 2D classification 2 rounds 2D classification 
(1,070,085 particles) (1,393,563 particles) 

3D classification 2 rounds 3D classification 


Mask on receptor Class 5, 89.9696, 246,392 particles 


9» Y 


Class 1 Class 2 Class 3 Class 4 Class 5 
9.28% 6.07% 62.62% 5.07% 17.02% 
88,533 P 59,130 P 698,935 P 47,740 P 175,747 P 


3D refinement (3.70 A) 


3D classification 
Without alignment 


Class 1 Class 2 Class 3 Class 4 Class 5 
4.27% 17.96% 5.93% 66.71% 5.14% 
17,946 P 74,311 P 24,884 P 283,561 P 21.437 P 
3D refinement (3.70 A) 
[ Transfer to Cryosparc 
Heterogeneous Refinement d 


GSFSC Resolution: 3.194 


102,581 particle 


Homogeneous Refinement 


Non-uniform Refinement 


(3.19 A) * 

558 È UNECTNN a ae ETT 2 

559 Supplementary Fig. 3 Cryo-EM data processing of CF102-A3AR-Gi complex. 

560 a. Representative image from cryo-EM dataset. Scale bar, 50 nm. b. Respresentative 2D 
561 average classification classes. Scale bar, 10 nm. c. Flow-chart of the cryo-EM data 

562 processing. d. FSC curves. e. The local resolution map. 

563 
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564 


565 TM1 TM2 TM3 TM4 TMS TM6 TM7 H8 a5 aN 
566 Supplementary Fig. 4 Respresentative regions of cryo-EM density maps are shown for 
567 the each transmembrane helical (TM) of AsAR and the a5 and aN helices of Gai, 
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570 Supplementary Fig. 5 Sequence alignment of adenosine receptors. The sequence 
571 alignment was generated with Jalview!4?! and depicts the N-/C-turminus, transmembrane 
512 helices (TMs), extracellular loops (ECLs), and intracellular loops (ICLs). Residues lining the 
573 orthosteric binding site are highlighted with red circles and annotated with GPCR 
574 Ballesteros-Weinstein numbering scheme. The C-terminus of the adenosine receptors were 
575 omitted. 
576 
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578 Supplementary Fig. 6 Effects of CF101 or CF102 on the A3AR mutants. These residues 
579 in AsAR formed hydrophobic interactions with the 3-iodophenyl group present in CF101 and 
580 CF102. 
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586 Supplementary Fig. 7 The orthosteric binding pockets among adenosine receptors. 
581 The positions highlighted indicate where unique residues occurred in ASAR compared to 
588 other adenosine receptors. The receptor names and their associated colors are shown 
589 above the models. The side chains in ASAR were depicted as bold sticks, while the 
590 corresponding side chains in other adenosine receptors were shown as thick sticks. 
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(7LD4) (6GDG) (8HDP) 
Supplementary Fig. 8 The binding cavities in adenosine receptors. 
a-e The binding cavities of the adenosine receptors are depicted as gray surfaces, with the 
bound ligands shown as sticks. The receptor names and associated PDB codes!*': ?^ 45] are 
indicated below each model. The unique subpocket in the A3AR as the red boxes circled. 
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597 Supplementary Table 1 Cryo-EM data collection, model refinement and validation 
598 statistics. 


A3AR-CF101-Gi complex A3AR-CF102-G; complex 


Data collection and processing 


Detector K3 K3 
Magnification 105,000 105,000 
Voltage (kV) 300 300 
Electron exposure (e-/À?) 50 50 
Defocus range (um) -1.0~-3.0 -1.0~-3.0 
Pixel size (A) 0.824 0.824 
Symmetry imposed C1 C1 
Initial particle projections (no.) 9,970,197 8,924,341 
Final particle projections (no.) 271,323 102,581 
Map resolution (À) 3.29 3.19 
Map resolution range (À) 2.60-4.40 2.40-4.20 
FSC threshold 0.143 0.143 


Model Refinement 
Refinement package 


PHENIX-1.17.1-3660 


PHENIX-1.17.1-3660 


Real or reciprocal space Real space Real space 

Model-Map CC (mask) 0.60 0.72 

Model resolution (À) 4.10 3.40 
FSC threshold 0.5 0.5 


B factors (À?, min/max/mean value) 
Protein residues 30.00/127.38/68.24 
Ligands 20.00/20.00/20.00 


30.00/135.93/68.27 
20.00/20.00/20.00 
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Model composition 


Non-hydrogen atoms 8,751 8,753 
Protein residues 1,126 1,127 
R.m.s. deviations 
Bond lengths (À) 0.001 0.005 
Bond angles (?) 0.398 0.900 
Validation 
MolProbity score 1.46 1.15 
Clashscore 8.51 3.62 
Rotamer outliers (96) 0.00 0.00 
Ramachandran plot 
Favored (%) 98.19 98.38 
Allowed (%) 1.81 1.62 
Disallowed (%) 0.00 0.00 


Data availability 
EMDB code 
PDB code 
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600 


601 Supplementary Table 2 Cell surface expression of A3AR and its mutants on CF101- and 
602 CF102-induced NanoBiT assay. 
ECso (NM)? pelleunace | 
expression (Relative 
CF101 CF102 to WT)? 
WT 3.74+0.7 1.50+0.4 100+2.5 
Y15!39A UD^ UD 113+5.0 
L90?32y 2.71+0.6 3.94+0.5 79.64+6.8 
L91?33A UD UD 107+12.1 
T94336A UD UD 91.5+4.7 
H95937A UD UD 121+7.1 
H95°37Q 7.7043.9 5.01+1.6 124+5.6° 
F168ECL2A >10000 UD 1164.0 
V169*C 2A NT* 1.07x0.3 85.9+6.0 
M174535A 1.973+0.5 NT 85.0+2.3 
M177938A »10000 UD 110+1.4 
S$181542A UD UD 1024+5.3 
S181542N UD UD 88.8+4.1 
1186547A UD UD 79.242.7 
11865477V 10.64+2.3 29.96*7.1 ^ 93.642.3 
W243°48A >10000 UD 65.5*6.2" 
L24695!A »10000 UD 81.0+4.9 
$247°5°H >10000 UD 82.5+6.2 
N250955A UD UD 78.2+5.0 
12536-58T 17.04+4.1” 13.11+4.5 90.1+9.8 
L264735A 4.08+0.8 5.53+1.5 73.743.8* 
1268739A UD UD 67.424.2" 
H272" ^A UD UD 12.74+3.6" 
603 
604 a Data shown are means + S.E.M. from at least three independent experiments. 
605 > UD indicates that the activation level is too low to determine EC»so values. 
606 * NT, not test. 
607 * P<0.01; **P<0.001 and ***P<0.0001 by one-way ANOVA followed by multiple comparisons 
608 test, compared with WT. 
609 
610 


611 
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612 Supplementary Table 3 Cell surface expression of A;AR/A2AAR/A?sAR and its relative 


613 mutant on CF101- and CF102-induced NanoBiT assay. 
CF101 CF102 Cell-surface 
max change max change expression 
Receptor i i : 
ECso(nM)? (Relative to ECso(nM) (Relative to (Relative to 
basal)? basal) WT)? 
AiAR 206.8+55 1.19+0.1 >10000 2.12+0.1 100+1.5 
A;AR-ECL3 161.9436 | 3.69:0.2" 176.64102  3.86+0.6" 91.14+3.2 
AzaAR 1280+727 1.53+0.1 UD» 1.03+0.1 100+2.9 
AzaAR-ECL3 74.29+16 2.21+0.1° 48.04+28 1.84+0.1 92.1+1.6 
A2BAR >10000 1.77+0.1 UD 0.95+0.1 100+7.6 
AzpAR-ECL3 1952+289 4.06+0.2~ >10000 2.72+0.2” 80.5+1.9" 
614 Se Oe hh c mh 
615 ? Data shown are means + S.E.M. from at least three independent experiments. 
616 > UD indicates that the activation level is too low to determine EC»so values. 
617 * P<0.01; **P<0.001 and ***P<0.0001 by one-way ANOVA followed by multiple comparisons 
618 test, compared with WT. 
619 
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